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ABSTRACT 

Characterization of the relatively poorly-understood progenitor systems of Type la supemovae is of great 
importance in astrophysics, particularly given the important cosmological role that these supernovae play. Ke- 
pler's Supernova Remnant, the result of a Type la supernova, shows evidence for an interaction with a dense 
circumstellar medium (CSM), suggesting a single-degenerate progenitor system. We present 7.5-38 /im in- 
frared (IR) spectra of the remnant, obtained with the Spitzer Space Telescope, dominated by emission from 
warm dust. Broad spectral features at 10 and 18 /im, consistent with various silicate particles, are seen through- 
out. These silicates were likely formed in the stellar outflow from the progenitor system during the AGB stage 
of evolution, and imply an oxygen-rich chemistry. In addition to silicate dust, a second component, possibly 
carbonaceous dust, is necessary to account for the short-wavelength IRS and IRAC data. This could imply a 
mixed chemistry in the atmosphere of the progenitor system. However, non-spherical metallic iron inclusions 
within silicate grains provide an alternative solution. Models of collisionally-heated dust emission from fast 
shocks (> 1000 km s" 1 ) propagating into the CSM can reproduce the majority of the emission associated with 
non-radiative filaments, where dust temperatures are ^80-100 K, but fail to account for the highest temper- 
atures detected, in excess of 150 K. We find that slower shocks (a few hundred km s" 1 ) into moderate density 
material (no ~ 50-250 cm" 3 ) are the only viable source of heating for this hottest dust. We confirm the finding 
of an overall density gradient, with densities in the north being an order of magnitude greater than those in the 
south. 

Subject headings: stars: AGB and post- AGB — ISM: supernova remnants — dust, extinction 



1. INTRODUCTION 

Kepler's supernova remnant (SNR) is the remains of the su- 
pernova (SN) of 1604 A.D., the third youngest known Galac- 
tic SNR, after Gl. 9+0.3 and Cas A. There is general consen- 
sus that Kepler resulted from a Type la SN, albeit one with 
an unusually dense circumstellar medium (CSM) (Blair et al. 
2007, and references therein; hereafter B07), evidenced by the 
extreme densities observed around Kepler (~ 2 — 100 cm" 3 ), 
which, at > 500 pc out of the galactic plane, should have sur- 
rounding densities of < 10 cm" 3 ( McK ee & Ostrikeri[l977h . 
as well as the nitrogen overabundance seen in optical spec- 
tra (IB lair et al.l 11991b . and high (~ 200 km s" 1 ) blueshift. 
Reynolds et al. (2007j) have suggested that Kepler may be a 
member of an emerging class of "prompt" Type la SNe, evolv- 
ing from the main-sequence to SN in < 5 x 10 8 y r. The 
distance to Kepler was determined by Sa nkrit et al.l (2005) 
to be 3.9 kpc, altho ugh other es timates range from 3 .3 kpc 
(Kat sudaet alj2008h to > 6 kpc dChiotellis et alJl20T3) . 

There is much debate in the general supernova literature 
about the progenitor systems of Type la SNe, particularly as 
to whether they result from single-degenerate systems (the 
explosion of a white dwarf that has accreted matter close 
to the Chandrasekhar limit of ~ 1 .4M from a companion) 
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or double-degenerate (the merger of two sub-Chandrasekhar 
mass white dwarfs). Ha emission from non-radiative shocks 
is seen at the periphery of the remnant, as we show in Fig- 
ure Q] In the case of Kepler, the high densities and presence 
of hydrogen in the CSM favor a single-degenerate scenario, 
since a double white dwarf system, even if it could eject a 
dense CSM, should not eject one rich in H. Furthermore, the 
presence of dust in the CSM implies that the companion star 
was an AGB star, since a main-sequence co mpanion should 
not pr oduce any significant quantity of dust. IChiotellis et alj 
d2012b suggested that Kepler resulted from a white dwarf ex- 
plosion in a symbiotic binary, and is expanding into the dense 
AGB wind from the progenitor system. 

As the only known Type la SNR encountering the wind 
from its pre-SN system, Kepler offers a unique opportunity 
to study the mass-loss from evolved progenitor stars. AGB 
stars are typically classified as O-rich or C-rich, with O-rich 
stars producing silicate dust and C-rich stars producing car - 
bonaceous dust in their outflows (Guh a Niyogi et all 1201 ll) . 
Through the interaction of Kepler's forward shock with this 
material, we can constrain the type of star present in the pro- 
genitor system. T he dust in Kepler is kn own to be dominated 
by silicate grains dDouvion et al.l 12001). implying formation 
in an O-rich AGB star. 

The formation of dust in O-rich stars is still not entirely un- 
derstood; here is one promising scenario. Large, transparent 
silicate grains form close to the stellar photosphere, within 
two stell ar radii. These are magnesium rich silicates, with- 
out iron (Norris et al. 2012). As these grains are driven out- 
ward by the stellar radiation pressure, temperatures fall and 
metallic iron condenses onto the grains, forming "dirty sili- 
cates" w ith substantial opaci ty at optical and near-IR wave- 
lengths ( Kemp er et al.l 12002b . This scenario has been con- 
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firmed in the laboratory by the identification of circumstellar 
grains among GEMS (glass with embedded metal and sul- 
fide s) grains; 1-6% of GEMS g rains are of circumstellar ori- 
gin dKeller & Messenge3l201 lh . However, the question then 
becomes one of how the presence of a white-dwarf companion 
can modify this process, both before and after the supernova. 

Prior to the explosion, accreting white dwarfs produce sig- 
nificant amounts of energetic radiation and shocks, quite pos- 
sibly enough t o dissociate molecules along the dust condensa- 
tion sequence (Buiarrabal et al. 2010). In an O-rich star, all C 
is locked up into CO, but this dissociation (either via radiation 
or shocks) can also affect CO molecules, leading to a mixed 
C/O chemistry in the outer dust formation zone. This has been 
observed in planetary nebulae (Guzman-Ra mirez et al.ll20~TTb . 
and carbonaceous grains have also been detected amongst 
silicate-rich dust i n the symbiotic binary system CH Cyg 
(ISchild et al.lll999h . 

After the SN, the forward shock from the explosion will 
race through the surrounding medium at several thousand km 
s , heating the gas to temperatures of > 10 7 K. This hot 
plasma will heat grains pres ent in the CSM, destro ying them 
via sputtering in the process dDwek & ArendtlH992l) . 

In B07, we reported Spitzer imaging of Kepler, confirming 
the presence of dense material ahead of the shock along with 
the existence of a density gradient in the north-south direc- 
tion, and placed an upper limit of < 0.1M Q of dust present. 
Here, we report followup observations with Spitzer 's Infrared 
Spectrograph (IRS), to study the processes listed above in the 
context of an SNR, and show how the data can be interpreted 
in the context of the wind from a pre-SN AGB star. Since 
Kepler's SNR represents the only known opportunity to study 
dust in the CSM surrounding a Type la SN, any characteriza- 
tion of the dust properties could ha ve an impact on the light 
echoes observed from Type la SNe (Patat et al. 20061). 

2. OBSERVATIONS 

We observed Kepler's SNR on 2009 March 29 with both or- 
ders of the long-wavelength (14-38 /im), low-resolution (LL) 
IRS module and order 1 (7.5-14 /im) of the short- wavelength, 
low-resolution (SL) instrument. Both have a wavelength- 
dependent spectral resolving power of A/AA = 64- 128. The 
entire remnant was mapped with the LL instrument, while se- 
lected regions were mapped with the SL module. The LL map 
consisted of four 30s cycles per position, with a step size in 
the perpendicular direction of 9" after each 120s observation. 
A total of 29 pointings in the perpendicular direction were 
necessary. We then shifted the slits in the parallel direction 
by 100" and repeated the process to cover the entire remnant 
with sufficient background. 

The three SL maps consisted of observations of varying 
depth to avoid saturation. The center-east knots were ob- 
served with 22 pointings, each with five cycles of 60s dura- 
tion, yielding a total map size of 56" x 43". The north rim po- 
sition was covered with 17 pointings (seven 60s cycles each), 
for a total map size of 56" x 33". The northwest position, 
by far Kepler's brightest region, was observed in 25 point- 
ings (five 14s cycles each), for a total map size of 56" x 49". 
A step size in the perpendicular direction of 1.85" was used. 
The outline of the SL maps is shown on Figure la. 

The LL map allows creation of a three-dimensional (two 
spatial and one wavelength ) data cube using the IRS con- 
tributed software CUBISm (S mith et al.ll2007l) . which allows 
spectral extraction from any rectangular region. The pixel size 
of the LL IRS instrument is 5.1", but Spitzer' s point-spread 



function is > 5.1" at wavelengths > 20 /im. Thus, we require 
a minimum region size of a 2x2 pixel square (10.2"xl0.2") 
when extracting LL spectra, with the additional requirement 
that any 2x2 pixel square falls within one IRS slit position, i.e. 
the two pixels do not fall on different slit pointings. For the 
background, we use an average of four off-source slit point- 
ings (two north of the remnant, two south). The mid-IR back- 
ground surrounding Kepler is highly uniform. 

Extracting a full spectrum from 14 to 38 /im requires 
matching of orders 1 & 2. The flux levels of these orders gen- 
erally did not match perfectly, despite being extracted from 
identical regions. However, since we are only interested in 
the shape of the spectra (see Section|3]l, we adjust the levels of 
the two orders to match, using the order overlap of 19.5-21.3 
/im. The relative offset, a product of the inherent calibration 
uncertainties for each order of IRS0, varied from 5% to 30%, 
depending on the region. 

The short wavelength end of LL order 1 (22-25 /im) con- 
tains a consistent noise pattern of regularly-spaced spikes, 
present in every spectrum. We used every Spitzer IRS data 
reduction program available, but all efforts to remove these 
fringes were unsuccessful. However, they have no discernible 
effect on the models fits discussed in Section [3] and we leave 
them in the data shown. 

3. RESULTS AND DISCUSSION 

The spectral variation from place-to-place in Kepler is sig- 
nificant. In Figure Q] we show both the 24 /im Multi-band 
Imaging Photometer for Spitzer (MIPS) image of Kepler and 
the Chandra X-ray image, overlaid with six different spec- 
tral extraction regions, the spectra from which are shown in 
Figure [2] The exact locations of these six regions are given 
in Table Q] Along the north and south rims (regions 3 and 
4), spectra show a dust-dominated continuum with no lines, 
while in the NW (region 2), lines from low-ionization state 
ions are present, the only place in the remnant where emis- 
sion lines are prominent. This corresponds to the brigh test re - 
gion of radiative shocks seen in the optical dBlair et al.lll991l) . 
Modeling of collisionally-heated dust grains in the postshock 
environment provides s trong constraints on the gas density 
(Bork owski et al.ll2006t IWilliams et al. 2006), but aside from 
a brief discussion in Section 13.21 we defer detailed spatial 
mapping of densities to a future publication. 

3.1. Dust Temperature 

Determining dust temperatures in SNRs is complex, since 
different grain sizes will be heated differently, with the small- 
est grains undergoing extreme temperature fluctuations on 
short timescales. However, we can make some simplifying 
assumptions and create a proxy representative of tempera- 
ture, in the same way that "color temperatures" are defined for 
stars. We created two narrowband images of Kepler in spec- 
tral regions away from lines or strong dust features (28 /im 
and 16.2 /im), and divided these images to make a ratio map. 
The value of the ratio at each image pixel was fit to a model 
of silicate dust at a single temperature, assuming a grain size 
of 0.05 /im. The output is shown in Figure lc. Grain "temper- 
atures" ranged from < 100 K to > 150 K, a high temperature 
for SNR dust. 

We made simple estimates of heating by UV photons from 
the radiative shocks propagating into modest CSM densities. 

6 http://irsa.ipac.caltech.edu/data/SPITZER/docs/dataanalysistools/tools/cubism/ 

7 http://irsa.ipac.caltech.edu/data/SPITZER/docs/irs/features/ 
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Specifically, we consider the case of a shock of 500 km s 1 en- 
countering a medium of «o = 250 cm -3 , and assume that 100% 
of the kinetic energy of the shock is converted into radiation 
(this clearly provides an upper limit, as such an approximation 
is unphysical). The radiation energy density from this shock 
would be ~ 1.2 x 10~ 9 ergs cm" 3 . Using the formulae pro- 
vided in lDraind (1201 lh . a radiation field of this magnitude will 
heat typical ISM grains of radius 0. 1 /xm to only ~ 50 K, far 
too low to account for dust seen in Kepler. Collisional heat- 
ing from hot gas in shock s is far more efficient in most SNRs 
(Draine & Salpeter 1979). There is a range of shock speeds 
in Kepler, from the < 100 km s" 1 radiative shocks in the NW 
to the > 3000 km s" 1 non-radiative shocks in the south rim. 
The latter, encountering the lowest density material («o ~ 1-5 
cm" 3 , B07) are capable of heating dust to 75-100 K. How- 
ever, if we assume pressure equilibrium, i.e., pv 2 is constant, 
then slower shocks of a several hundred km s" 1 should be en- 
countering material with densities of 50 to se veral x 100 cm" 3 . 
We use the grain heating models de scribed inlBorkowsk i et all 
(2006) and IWilliams et al.l d2011l) to examine the effects of 
collisional heating on grains in these slower shocks. Briefly, 
grains are heated via collisions with hot ions and electrons in 
the post-shock gas, with final grain temperatures being a func- 
tion of both the plasma temperature and density. Shocks such 
as these can collisionally heat grains to > 150 K. 

3.2. North-South Density Gradient 

Kepler's north limb is significantly brighter than the south 
across the EM spectrum. In B07, we attributed this to the 
northern shock e ncountering material 4 -9 times denser than 
that in the south. Katsud a et alJ ([2008) measured the proper 
motions of X-ray emitting filaments around the periphery, 
finding that shock velocities in the northern filaments were 
1 .5-3 times slower than in the south, implying a density ratio 
of 2-9. 

We extracted IR spectra from near the locations where 
proper motions were measured. Spectra of collisionally- 
heated dust are most sensitive to the post-shock gas density, 
and by making reasonable assumptions about the ion and elec- 
tron temperature, we can fit a spectral model where density is 
the only free parameter. We assume pure silicate dust and 
use a x 2 minimization algorit hm to fit the spectra from 2 1 to 
33 /im, which we have found (Will iams et al.ll201ll) to be the 
most reliable region of the spectrum to fit to obtain a gas den- 
sity. In the non-radiative sections of the north rim (region 3), 
we use the values derived in B07 of T, ■ = 8.9 keV and T e = 1 .4 
keV, obtaining a post-shock density of «h = 42 (41,43.5) 
cm" 3 , with a reduced \ 2 of 1.3 (values in parentheses are the 
90% confidence limits, but represent only the statistical errors 
on the fit). In this model, ~ 60% of the dust is destroyed via 
sputtering, roughly consistent with our initial modeling of the 
MIPS data in B07. 

In the southernmost rim (region 4), X-rays are dominated 
by nonthermal emission, and determining the plasma tem- 
perature is difficult. Using the same temperatures as above, 
we find «h = 4 cm" 3 . However, it is likely, given the high 
shock speeds and low densities, that far less temperature 
equili bration has taken place between the ions and elec- 
trons (Ghavami an et al.l 12007). If we thus assume an ion 
temperature of 50 keV, allowing for heating of electrons 
through Coulomb collisions alone, we derive a density of 
n H = 2.0 (1.8,2.2) cm" 3 (reduced X 2 / = 0.99). Only 10% of 
the dust is destroyed via sputtering in this region, owing to 



the much lower gas densities. We show the fits from both re- 
gions, 3 & 4, in Figure [3] Slightly farther south, in the region 
that in X-rays is dominated b y thin, synchrotron emitting fila- 
ments (Reynold s et alJl2007l) . we find a weak dust continuum 
that is best fit with a density of 1 cm" 3 . At Spitzer's reso- 
lution, it is difficult to determine if this continuum is real or 
simply due to the wings of the point-spread function from the 
brighter emission just above. Thus, we consider n# = 1 cm" 3 
to be an upper limit for the southernmost rim. The addition of 
graphite grains to these models (see Section |4) increases the 
derived densities by ~ 20%, but the overall density gradient 
is unchanged. 

3.3. Silicate Dust 

Strong features, most easily attributable to silicate dust, 
exist in the spectra shown in Figures 2-5 at ~ 10 and 18 
/im. In Figure @] we show the total 7.5-40 /im spectrum 
from region 3, overlaid with two different model fits to the 
data: a simple blackbody model (T ~ 120 K) and an emis- 
sion model for graphite grains. The broad emission fea- 
tures, visible in the data at ~ 10 /im and ~ 18/xm are eas- 
ily seen in excess of either of these two models. While 
these silicate features are often observed in absorption in cold 
dust of the ISM, they are routinely observe d in emission in 
places where dus t is warm, such as SNRs (Rho et al.l 120081; 
iDwek et ail 1201 Oh or the optically-thin d ust shells of O-rich 
AGB s tars with moderate ma ss-loss rates dGuha Niyogi et all 
1201 ll) . iDouvion et al.l (120011) first identified these features in 
Kepler in Infrared Space Observatory spectra, reporting that 
they could be fit with stan d ard si licate dust with optical con- 
stants from iDraine & Led (119841) . We find no evidence for 
silicate absorption features anywhere in the remnant. 

On the smallest measurable scales (10.2" x 10.2"), a pecu- 
liar feature is present in the spectra of about two dozen regions 
within the remnant. In these regions, the 18 /im silicate fea- 
ture appears to be so strong that it peaks above the level of 
the continuum, dipping back down sharply on the long wave- 
length side. An example of a spectrum like this is shown in 
the upper left panel of Figure [2] We have no complete ex- 
planation for this feature. When slightly larger spatial extrac- 
tion regions are used, the feature disappears. We examined 
the individual spectral observations that went into making the 
cube using SPICE, the Spitzer Custom Extraction Software, to 
search for evidence of the feature's existence, and found that 
the appearance and disappearance of the feature depended en- 
tirely on the spectral extraction method used. The effect is 
strongest in regions of the remnant where strong gradients in 
the IR surface brightness exist. Because of this, we believe 
it is most likely that it is an artifact of the spectral extraction 
algorithms used by CUBISM on the smallest spatial scales, 
and not an absorption feature or abnormally strong unknown 
emission feature. 

3.3.1. Possible Ejecta Dust? 

Two possibilities exist for the dust seen in Kepler: newly 
formed ejecta dust or dust created in the AGB wind of the pro- 
genitor system. Ejecta dust is typically associa ted with core- 
collap se SNe, but recent theoretical work by iNozawa et al.l 
(2011) has shown that it may be possible for up to 0.2 M Q 
of dust to form in the ejecta of Type la SNe. However, to 
date, there have been no d etections of ejecta d ust firmly asso- 
ciated with Type la SNe. iGomez et al.l d2012l) report a small 
amount (~3 x 10~ 3 M®) of dust detected in Herschel obser- 
vations (comparable to what we found in B07), and mention 
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the possibility that they cannot rule out, based on their Her- 
schel data alone, a contribution from ejecta dust, but they do 
not claim detection. 

Since the first detection of ejecta dust in a Type la SNR 
would be of extreme interest, we have searched extensively 
for any signs o f it in Kepler. We h ave gone about this in 
two ways. First, Burkev et al. (2012) separate the X-ray data 
from Chandra (Rey nolds et al.ll2007h into regions dominated 
by ejecta emission, thermal CSM emission, and nonthermal 
synchrotron radiation. A detailed morphological comparison 
between these data and ours shows no obvious evidence for a 
correlation between IR emission and X-ray ejecta dominated 
regions. While we cannot rule out an ejecta dust contribu- 
tion to the IR spectra (complications from projection effects 
cannot be ignored), it does not appear to be significant. 

Secondly, another comparison can be made between the lo- 
cation of dust emission and the emission line maps produced 
by CUBI SM, a technique used for a spectroscopic study of 
Cas A bv lRho etaH d2008l) . where they use [Fe II] and other 
emission lines. Kepler is ~ 20 times fainter than Cas A and 
the maps produced are of much lower signal-to-noise, and the 
[Fe II] lines, which can come from either ejecta or shocked 
or photo-ionized gas are the only ones strong enough to make 
a spectral map. The [Fe II] emission line maps at 17.9 and 
26 /im show clear emission from several regions in Kepler, 
but again, no obvious correlation exists between the regions 
of strong line emission and peculiar dust spectral features, as 
is the case in Cas A. Again, while we cannot rule out an ejecta 
dust component to the IR spectra, we unfortunately have no 
direct evidence for it. 

Lastly, as we point out in B07, there is a clear morphologi- 
cal similarity between dust emission seen with all of Spitzer's 
instruments (the Infrared Array Camera; IRAC, MIPS, and 
IRS) and the Ha emission shown in Figure [TJ Since this Ha 
emission is associated with the non-radiative shocks at the 
shock front, we conclude that the dust is associated with the 
swept-up CSM. 

3.3.2. CSM Dust 

It thus appears that the shock is primarily encountering dust 
formed in the outflow from the progenitor system, either the 
precursor star itself or the donor. AGB stars produce a slow, 
dense wind that is ideal for the formation of grains. The ther- 
modynamic condensation sequence for dust formation in the 
atmosphere of an O-rich AGB star is complex and not fully 
understood (see Section [TJ, but the end result in high mass- 
loss rate outflows is a significant a mount of amo rphous sili- 
cate dust that we also see in Kepler (Waters 201 1). 

While it is clear that projection effects ensure that the spec- 
tra in Kepler sample regions of differing temperature, density, 
and composition, we can nonetheless arrive at some general 
conclusions even without detailed modeling. First, the spec- 
tra from the central knots, shown in Figure |2] are remarkably 
similar to the spectrum from the north rim, consistent with the 
hypothesis that emission from these central knots results from 
CSM dust. Second, although the overall shape of the contin- 
uum at long wavelengths is relatively well-fit, the 18 /im sili- 
ca te feature is no t well-f it with the astronomical silicate model 
of Draine & Led (119841) . which places the centroid of the fea- 
ture at ~ 18.2 /im. The silicate features seen in the spectra of 
Kepler in the north, northwest, and central knots are centered 
at ~ 17 /im. This is most easily seen in Figure|3] where a fit to 
the long-wavelength spectrum of region 3 significantly under- 
predicts emission from 15-18 /mi. To further test this conclu- 



sion, we follow the technique of Guha Nivog fet al.l d201 lb of 
fitting a blackbody model to the underlying continuum (using 
a least-squares routine) and dividing this out of the spectrum. 
What remains is the absorption efficiency of the silicate dust 
grain as a function of wavelength. While the 10 /im feature 
is in the correct place (to the extent that our short-wavelength 
data will allow us to determine this), the 18 /im feature is not, 
and is much closer to 17 /im. We find no dust tem perature 
at which the optical constants of iDraine & Led (1 19841) can ac- 
count for the features observed. This may point to a different 
chemistry within the silicate grains. This is perhaps not en- 
tirely surprising, as there are sy stematic differences be tween 
CSM and ISM silicates features dOssenkopf et al.ll 19921) . Fur- 
thermore, a great variatio n in silicate features is seen in AGB 
outflows (Henning 2010), suggesting marked differences in 
physical and chemical properties among CSM grains. We 
also note that a similar mismatch of the 18 /nm feature wa s 
observed in the central region of M81 (SmithetaLllIOlO), 
where the feature peaks at 17.2 /im. The authors there found 
a better (but far from perfect) match when using the dielectric 
fu nction of amorphous o livine as measured in a laboratory 
by lDorschner et all (fl 995). But distinguishing among various 
types of silicates in Kepler is a difficult task that requires sep- 
arate efforts involving more sophisticated spectral modeling 
and reliance on mo re extensive laboratory measurements (see 
ISpeck etakl (1201 ll) for recent experimental studies of amor- 
phous silicates). 

We find no indication of the strong spectral features pro- 
duced by crystalline silicate particles, such as those seen in 
the spectra aro und cool, evolved O-rich stars (Morris! 120081 
Henning 2010). These features are generally most prominent 
in high mass-loss rate (> 1O~ 5 M0 yr" 1 ) AGB outflows, so it 
is possible that the mass-loss rate of the progenitor system 
of Kepler was not extreme. But crystalline silicates may be 
common as well at low mass-loss rates, just harder to detect 
(e.g., see de Vries et al. 2010 and references therein), and the 
signal-to-noise ratio in Kepler may make any such detection 
impossible. Crystalline silicates are rarely detect ed in SNR 
for a c ounter-example where they are detected, see iKoo et all 
201 ll) and compare the spectra from MSH 15-52 to that from 
Kepler), and a likely explanation for this is that impacts with 
thermal ions amorphize crystalline silicates, changing their 
chemical structure. For the fast non-radiative shock in the 
north of Kepler, with 7} = 8.9 keV and ~ 40 cm" 3 , the 
time required for the complete amorphization is only several 
years for grains ^0.1 /i m in radius, based o n critical proton 
and a particle doses in iGlauser et al.l (120091) . In the south, 
densities are an order of magnitude lower than in the north, 
and the amorphization timescale is inversely proportional to 
density. Crystalline silicates might survive under these con- 
ditions, but again, would be harder to detect because of low 
IR fluxes. Crystalline material might also survive in cores of 
large grains, shielded from energetic ions by overlying amor- 
phized surface layers, particularly in slow shocks driven into 
the dense CSM where ions do not hav e enough energy to pen- 
etrate into grain cores. Glaus er et al.l (1200 9) showed that such 
partially crystalline grains can still produce detectable crys- 
talline emission features, but we see no evidence for them in 
Kepler. 

4. THE NEAR-IR EXCESS 

Since dust in Kepler exhibits prominent silicate emission 
features, it must have condensed in a stellar outflow with dom- 
inant O-rich chemistry (C/O ratio < 1). However, perhaps 
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the most intriguing spectral feature of Kepler is that while the 
spectra are clearly dominated by silicate grains, these alone 
cannot account for the short-wavelength data, including the 
IRAC fluxes. We use IRAC data obtained as part of our ear- 
lier survey of Kepler with both IRAC and MIPS (PI W. Blair, 
program ID 3413, see B07). We show this in Figure [5] the 
IRS short-low spectrum from the bright NE knot (rectangu- 
lar region in Figure Id). This NE knot is not unique; other 
locations in Kepler show a similar near-IR excess. A local 
background, subtracted from both sides of the knot, was used 
to highlight the hottest dust. For the IRAC data, we measured 
the flux from an identical region (CUBISM allows creation of 
region files that can be read by ds9) at 5.6 and 8.0 /im. The 
8 /im flux was roughly consistent with the IRS data, and thus 
does not provide an additional constraint, but the 5.6 /im flux 
provides a strong confirmation of the near-IR excess. At 3.6 
and 4.5 /im, Kepler is too faint and the knot is undetected. 

A model fit to the data with pure silicates underpredicts the 
IRAC 5.6 /im emission by two orders of magnitude. This 
near-IR excess cannot be explained by synchrotron emission, 
nor by an additional hot silicate component. This near-IR ex- 
cess has been a long-standing issue in observ ations of both the 
ISM and envelopes around late-type stars (Ossenkopf et al. 
1992). An additional component of graphite dust provides 
a better explanation than pure silicates, as shown in Fig- 
ure |5] (though the fit is far from perfect at the shortest wave- 
lengths), which results from a shock model (v, = 500 km s -1 , 
«h = 250 cm" 3 ) of collisionally-heated dust, where silicates 
and graphites are assumed to have a power-law size distribu- 
tion (a = -3.5) from 1 nm to 0.5 /im and a mass ratio of sili- 
cates to graphites of two. Such a model, which still underpre- 
dicts the observed emission, requires a mixed C/O chemistry 
in the AGB progenitor system . This is possible in rar e S-type 
AGB stars where C/O rs 1 dFerrarotti & Gaill 120021) . Mixed 
C/O chemistry can operate in S stars if shocks are generated 
within their winds by stellar pulsations, le ading to formation 
of both carbonaceous and silicate grains dCherchneffl 120061 
Smol ders et al.ll20Tol) . If dust in Kepler was produced by an 
S-type AGB star, this might explain the presence of silicates 
together with carbonaceous grains. Unfortunately, we are un- 
able to determine the C/O ratio from UV and X-ray observa- 
tions because of the high interstellar absorption towards the 
remnant. 

We also consider more common M-type (C/O < 1) AGB 
outflows, where in single stars nearly all C is locked in CO 
molecules. In a single-degenerate Type la progenitor sys- 
tem, an accreting white dwarf produces energetic UV/X-ray 
radiation, and thermonuclear outbursts drive shocks into the 
ambient CSM. A fraction of CO molecules can be dissoci- 
ated by shocks or energ etic photons in the dust -forming re- 
gions in these systems (Bujarrab al et al.l feOlOl) . This cre- 
ates favorable conditions for formation of carbonaceous dust 
even if C/O < 1. While quantitative models of carbona- 
ceou s dust formation in O-rich sym biotic systems are lack- 
ing, iGuzman-Ramfrez et al.l (1201 ll) demonstrate how ener- 
getic photons can lead to mixed carbon-oxygen chemistry in 
O-rich outflows of Galactic bulge planetary nebulae. Their 
models show that formation of hydrocarbon molecules, pre- 
cursors to PAHs and hydrogenated carbonaceous dust, is more 
efficient at larger C/O ratios. 

However, the near-IR excess could also be produced by 
metallic iron (or various iron oxides). A similar near-IR ex- 
cess, with respect to models of pure silicate dust, was found 



bv lKemper et al.l d2002l) in their study of the O-rich AGB star 
OH 127.8+0.0. They examined the chemical abundances of 
various materials likely to be present in AGB atmospheres, 
considering several potential (non-carbonaceous) dust species 
that could form in addition to silicates, finding that non- 
spherical metallic Fe inclusions within the silicate particles 
of ~ 4% (by mass) could best account for the 3-8 /im con- 
tinuum. Such grains have recently been identified in terres- 
trial laboratory samples by Keller & Messenger! (1201 ll) . who 
conclude, based on their isotopic compositions, that some of 
these grains could originate from A GB stars. By compa ring 
the dust mass absorption coefficients dKemper et al.ll2002h for 
these "dirty silicates" to the near-IR excess we see in Kepler, 
we find that a similar amount of non-spherical Fe-inclusions 
(~ 5% by mass) in amorphous silicate particles can account 
for the IRAC 5.6 /im flux. This model does not require a 
mixed C/O chemistry in the progenitor system. This is also 
true of other possible models with separate po pulations of 
non-si licate Fe-bearing grains (e.g., see models of lDwek et al.l 
(2010) for SN 1987A which shows a similar near-IR excess). 

5. SUMMARY 

While other remnants of Type l a SNe at comparable ages 
to Kepler are interacting with ISM ( Williams et al.l2011l) . Ke- 
pler is unique in interacting with CSM, providing a laboratory 
to study the interaction of the forward shock with the mate- 
rial ejected by the pr ogenitor system during its AGB phase 
(Chiot ellis et al J 120121) . While Kepler is roughly spherical in 
shape, we find more than an order of magnitude density con- 
trast between the northern and southern rims. The strong sil- 
icate features observed here are indicative of the wind from 
an O-rich AGB star. However, there is clear evidence in 
the short-wavelength emission for an additional dust compo- 
nent, with either carbonaceous dust or metallic-Fe inclusions 
within amorphous silicates providing a possible explanation. 
The presence of carbonaceous grains would imply a mixed 
oxygen-carbon chemistry in the progenitor system. There 
is no morphological evidence for newly-formed ejecta dust. 
Small scale studies of the silicate features seen here, the 9.7 
/im feature in particular, will be an excellent observational 
goal for future IR telescopes, such as the James Webb Space 
Telescope, which can search for crystalline silicates and help 
distinguish between carbonaceous and iron grains. 
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TABLE 1 

Spatial Locations of Spectral Extractions 



Region R.A. (J2000) Dec. (J2000) Size (arcsec) 



1 17:30:43.1 -21:29:30.9 10.2 x 10.2 

2 17:30:36.0 -21:28:45.5 20.4 x 20.4 

3 17:30:44.1 -21:28:00.7 25.5 X 20.4 

4 17:30:40.4 -21:31:06.2 56 x 20.4 

5 17:30:39.9 -21:29:24.1 20.4 x 20.4 

6 17:30:42.7 -21:29:46.6 30.6 x 20.4 



NOTE. — Positions listed are the center of the region. Sizes of regions are given as length x height. 




109 113 116 120 123 127 130 134 138 

FIG. 1 . — a) Spitzer 24 fim image of Kepler (B07), o verlaid with re gions corresponding to spectra shown in Figure 2. Dashed white boxes represent locations 
of IRS SL maps, b) Three-color Chandra X-ray image ( Reynolds et al. 2007) with identical spectral regions overlaid, c) Temperature map (in K), as described in 
text, with scale indicated by color bar at the bottom, d) HST/ACS image with Ha (red), [N II] (green), and [O III] (blue). Regions appearing white are bright in 
all bands, indicating radiative shocks. Cyan box shows extraction region for spectrum in Figure 5, discussed in section 4. An outermost contour from just outside 
the forward shock is shown in panels c and d. N is up and E to the left in all images. Kepler is ~ 3.8' in diameter. 




FIG. 2. — IRS long-wavelength spectra of selected regions in Kepler, chosen to highlight the variation in dust throughout the remnant. Spectra are color-coded 
to regions shown in Figure 1 and plotted on a linear scale in surface brightness, averaged over the region. 
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15 20 25 30 35 

A (//m) 

FIG. 3. — LL (~ 14-38 fim) spectra from regions 3 (green) and 4 (yellow), overlaid with dust model fits as described in the text. Region 4 is multiplied by a 
factor of 40 to be shown on the same scale. The fits assume pure silicate dust and are made only to the 21-33 (im region of the spectrum. From the fit to region 3, 
it ca n also be seen that the 18 /im silicate feature cannot be adequately reproduced by these models, which assume the optical constants of "astronomical silicate" 
from lDraine~& Lee ( 198|). 
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FIG. 4. — The SL & LL (7.5-40 ^m) spectrum of region 3, plotted on a logarithmic scale. A single-temperature blackbody model and a single-temperature 0.05 
/mi graphite grain model are fit to the data and overlaid to show the existence of the broad 10 and 18 (im silicate emission features. 
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FIG. 5. — SL (7.5-14 fim) spectrum (in surface brightness, on logarithmic scale) of NE emission knot (cyan region in Figure Id), overlaid with shock model fit 
as described in text. The IRAC 5.6 fim data point, used as a constraint to the fit, is shown as a star. 



